Incorporating fermions in the Gaussian Effective Potential: the Higgs-Top Sector 
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The Higgs-Top model is studied by a non-perturbative variational extension of the Gaussian 
Effective Potential that incorporates fermions. In the limit of a very strong Yukawa coupling the 
one-loop result is shown to follow a single-parameter scaling while the gaussian fluctuations give 
rise to important deviations from scaling and to a reduction of the vev and of the top mass. A 
good general agreement is found with lattice data when a comparison can be made. The vacuum is 
shown to be stable for any choice of the Yukawa coupling, at variance with renormalized perturbation 
theory. Analytical results are provided for few observables like the renormalized mass of the Higgs 
boson and its wave function renormalization constant. Extensions to gauge theories like QCD are 
briefly discussed. 
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I. INTRODUCTION 

Above the QCD scale, when all the couplings are small 
enough, the Standard Model (SM) of the fundamental in- 
teractions can be studied by perturbation theory. How- 
ever, there are energy ranges and sectors of the SM where 
non-perturbative methods are welcome for comparison 
and control of the perturbative approximation. For in- 
stance when the QCD becomes strongly coupled, or in 
the hypothesis of a large Higgs mass and a strong self 
coupling of the Higgs sector, even if almost ruled out by 
experiments. An other interesting system is the Higgs- 
Top sector, where non-perturbative effects could be found 
because of the Yukawa interaction y sa 0.7 which is not 
very small as for other fermions. 

The most widely used non-perturbative approach is 
the numerical simulation on a finite lattice, which allows 
for an exact treatment of the Lagrangian, but is plagued 
by many shortcomings like the small size of the sample, 
border effects, the lack of any analytical result and some 
specific problems about incorporating fermions [3-0|- F° r 
instance a duplication problem does not allow a simple 
description of a model with an odd number of fermions 
or even one only fermion[7j. Moreover the finite lattice 
spacing is equivalent to the existence of a finite energy 
cut-off that cannot be removed. 

This last point does not seem to be an important is- 
sue any more, since our modern understanding is that 
the SM must be considered as an effective model holding 
up to a finite energy scale. The triviality of the scalar 
theory requires that a finite cut-off should be present in 
the model, and even a large cut-off can be simulated on 
the lattice if the sample size is large enough and allows 
for long correlation lengths compared to the lattice spac- 
ing. On the other hand, the existence of a finite cut-off 
has opened the way to a class of variational approxima- 
tions, since their typical UV problems are cured by the 
existence of a limited energy range. The Gaussian Ef- 
fective Potential (GEP) [si— Fill] is a variational tool that 
has been recently shown to provide a very good agree- 



ment with lattice data whenever a finite cut-off is used 
in the analytical derivation[T2|. The GEP is not exact, 
but is based on a non-perturbative approximation, and 
its predictions are expected to hold even when pertur- 
bation theory breaks down. Moreover the GEP yields 
simple analytical results and has been recently shown to 
sheld some light on the physics of several systems r ang ing 
from superconductors (l3ir3 to magnetic systems |l5j| to 
non-Abelia n gaug e theories [16[ and the Higgs sector of 

the SMpl [rrni9|. 

It would be interesting to extend the GEP to the Higgs- 
Top sector of the SM, which would require the inclusion 
of fermions in the derivation of the effective potential. 
The problem is of some interest by itself because of the 
failure of incorporating fermions that was reported in the 
past. In fact a direct attempt to include fermions was 
shown to give a result that is equivalent to the perturba- 
tive one-loop fermionic term of the effective potential |20j. 
Thus for fermions the GEP has always been regarded as 
useless. Quite recently an hybrid method has been pro- 
posed for incorporating fermions in the GEP[15]: the 
technique has been tested in the framework of the well 
studied two-dimensional Hubbard Hamiltonian at half 
filling, describing a gas of interacting fermions in the 
antiferromagnetic phase. The method predicts the ex- 
act magnetization in the strong coupling limit, and for 
large couplings it improves over other approximations like 
RPA, in close agreement with lattice data. Instead of at- 
tempting to evaluate the GEP directly, in this method 
the fermions are integrated exactly, and the resulting ef- 
fective Lagrangian is expanded in powers of the scalar 
field. Then the GEP is evaluated by the usual varia- 
tional method. In the strong coupling limit the second 
order term of the expansion has been shown to be enough 
for an accurate description of the fermionic fluctuations 
that are included at the Gaussian level(l5|. 

In this paper the same hybrid technique is used for in- 
corporating fermions in a self-consistent variational ap- 
proach to the Higgs-Top model, a toy model containing a 
self-interacting scalar particle interacting with a massless 
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fermion through a standard Yukawa interaction. In this 
model, that mimics the mechanism of mass generation of 
the SM, the symmetry breaking is driven by the Yuakawa 
interaction, and even when the classical scalar potential 
does not show any symmetry breaking, the effective po- 
tential can have a minimum with a non vanishing vacuum 
expectation value (vev) of the scalar field. Thus while 
the fermions acquire a mass through the Yukawa interac- 
tion, it is the fermion that drives the symmetry breaking 
determining the vev. A naive discussion of symmetry 
breaking, based on the classical potential, does not work 
in the present model where the vev must be determined 
by the full self-consistent quantum effective potential. 

In the next sections the GEP is evaluated for the Higgs- 
Top model, and the result is compared with available lat- 
tice simulation data[21] and with the standard one-loop 
approximation. A full agreement is found with the lat- 
tice data when a comparison can be made. Moreover 
the method is almost analytical, and explicit analytical 
results can be obtained for the pole of the scalar prop- 
agator, the wave function renormalization constant and 
other relevant observables. For instance, by inspection 
of the analytical expression of the effective potential, the 
vacuum can be shown to be stable for any choice of the 
bare parameters, at variance with renormalized perturba- 
tion theory. The method is non-perturbative, and devi- 
ations from the one-loop approximation can be regarded 
as a measure of non-perturbative effects that are shown 
to be large in the strong coupling limit. 

Besides the physical relevance of the present model, the 
problem of incorporating fermions in a non-perturbative 
method is by itself relevant because of the possible ex- 
tension to non-Abelian strongly interacting theories like 
QCD. For instance the GEP has been evaluated for the 
pure SU(2) theory [16| but no previous attempts of in- 
corporating fermions have been reported . The present 
method could be explored in that context. 

The paper is organized as follows: in Section II the 
method is described in detail; in Section III the GEP 
is compared with lattice data; in Section IV the strong 
coupling limit is explored and compared with the 1-Loop 
result; some final remarks and directions for future work 
are reported in Section V. 



II. GEP IN THE HIGGS-TOP MODEL 

In the Euclidean formalism the Higgs-Top model is de- 
scribed by the Lagrangian 

£> = £</> + £>t + £>y (1) 

where is the Lagrangian of a self- interacting scalar 
field <j) 

-c* = ia / V0 /1 + TO) ( 2 ) 



with a classical potential 

VM = \™ 2 B 4> 2 + ^A B 4 , (3) 

Lt is the Lagrangian for a set of Nf massless free Fermi 
fields ipj 

N f 

L t = - i J2^1 >1 d^ J (4) 

3=1 

and L y contains a set of Yukawa couplings 

L y = S )Zy^^i$- ( 5 ) 

3=1 

Here we take yj — y so that the set of Fermi fields ipj 
can be regarded as a gauge multiplet, and for Nf = 3 the 
model describes a quark interacting with a scalar field. 
If the symmetry is broken and the scalar field has a non- 
vanishing vev, v =< (p >, then the Higgs field h is defined 
by the shift h = <p — v and the fermion acquires the mass 
to = yv. We refer to m as the top mass even if the present 
study applies to any fermion with a Yukawa coupling to 
the scalar field. 

Let us follow the method of Ref.[l5j and integrate the 
fermions exactly. The exact action can be written in 
terms of a shifted Higgs field h = <p — <p where <p is a 
generic constant shift. The field h becomes the standard 
Higgs field when ip = v, but we leave the variable <p 
unconstrained at this stage. Integrating out the fermions, 
the bilinear fermionic terms are replaced by an effective 
action 



S = J Ld 4 x^ J L h d 4 x + S eff [h] (6) 

where 

£ h = ^hd^h + V c (<p + h) (7) 
and the effective action S e ff is given by 

S eff [h] = -N f log dct [G- 1 + yh] m=w (8) 

with the fermionic inverse propagator G" 1 defined as 

G-^-i^ + m (9) 

and the mass set to the value m — y(p. These equations 
are exact and hold for any choice of the variable ip as 
a consequence of the exact integration of the fermionic 
fields. 

It is obvious that the effective action in Eq.© is far 
too complicated to be treated exactly. From a formal 
point of view it can be expanded in powers of the field h 
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S eff [h] = -N f TrlogG- 



-N 



i 



l y n Tr{(G m -h) n } 



(10) 

where the products G rn ■ h ■ G m ■ h ■ ■ ■ are non local prod- 
ucts of operators with the internal space variables inte- 
grated over. 

The expansion is convergent if the norm ||y/iG m || < 1. 
We will show that in the strong coupling limit y 3> Xb ~ 
1, in units of the cut-off (or of the lattice spacing), the 
physical range of parameters, where m <C 1 is reached 
for 0.2 < ms/y < 0.25. In fact in that limit the one- 
loop approximation is entirely characterized by the sin- 
gle parameter ms/y which takes its critical value when 
m B /y = i/N]/(2tt), that is w 0.28 for N f = 3. We will 
discuss this scaling later, in Section IV. In that physical 
range, provided that we do not reach the critical point, 
the average < hh Z/M R where Z « 1 is a wave func- 
tion renormalization factor and Mr is a renormalized 
mass of order Mr k, tub, while in units of the cut-off 
< G m G m >« -l/(47r 2 ), so that 



\\yhG n 



y 



Tr < hG m hG m >f 



( yz 



(11) 



and for Z = 1, Mr, w tub we see that the expansion 
makes sense provided that 



\yhG r , 



0.08 • — < 1 

m B 



(12) 



In the physical range of parameters where ms/y 0.25 
this is never too large, even in the very strong coupling 
limit. 

In fact the expansion in Eq. fTU]) is not a perturbative 
expansion in the parameter y, but can be regarded as 
an expansion in terms of the fluctuating field h. Pro- 
vided that < hh > is small, the gaussian fluctuations 
are enough and the quartic and higher order terms can 
be dropped without affecting the non perturbative na- 
ture of the approximation. This is the case for the Hub- 
bard model of antiferromagnetism in the strong coupling 
limit [l5|]. The approximation breaks down at the critical 
point where the fluctuations of the field h are large and 
< hh > diverges. However the triviality of the scalar 
theory requires that a finite cut-off is retained, that is 
equivalent to say that the critical point is never reached. 
In the physical range of the parameters, before the crit- 
ical point is reached, even in the strong coupling limit 
the fluctuations are small and we can neglect the higher 
order terms in the expansion. This approximation spoils 
the variational character of the method, but the non- 
perturbative nature of the approximation is maintained 
as for RPA approximation that is very similar to the 
present method from a formal point of view. In principle 
the approximation can be improved by inclusion of higher 
order terms in the expansion, yielding a more complicate 
set of coupled equations that can be solved by numerical 



techniques. However the first non- vanishing contribution 
comes from the quartic term < hhhh > that is much 
smaller than the second-order term in the physical range 
of the parameters. Its actual value can be evaluated for 
a more accurate control of the approximation. 
The truncated expansion for S e ff reads 



Sj$ f [h] = -NfTrlogG^ - yN f Tr(G m ■ h)- 
+ ^y 2 N f Tr(G m ■ h ■ G m ■ h) 



(13) 



and our starting point is the vacuum-to-vacuum ampli- 
tude 



Z = D h e 



C h d*x + S<$ f [h] 



(14) 



Inserting a gaussian trial functional, the amplitude Z is 
written as 



Z= \D h e 



h(x)g 1 (x 1 y)h(y)d 4 xd 4 y + S int [h} 



where g 1 (x,y) is a trial inverse propagator, and 



Sint[h] = / h(x)g 1 (x,y)h(y)d 4 xd 4 y+ 



(15) 



■J L h d 4 x + S { e f f [h} 



(16) 



that is regarded as an interaction term. We define the 
gaussian average of a generic operator according to 



1 



(o) = y I Vh0e 

where 



Zo = / V h e 



h{x)g 1 (x,y)h(y)d 4 xd 4 y 



h(x)g 1 (x,y)h(y)d 4 xd 4 y 



(17) 



(18) 



It follows that < h >= and < h(x)h(y) >= g(x,y). 
Moreover < cf> >— <p and the variable <p is the expectation 
value of the scalar field. 



By Jensen inequality |16l . |22j the effective potential 

(19) 



e qua 

V(<p) is bounded[llJ, [23J by a gaussian functional 



V(tp) < V GEP [g] 
where the gaussian functional is defined as 
VcEp[g] = (Sint)+ 

K^g' 1 ^, y)Hy)d 4 xd A y 



log D h e 2 



(20) 
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and is equivalent to the first-order effective potential in 
presence of the interaction term Si nt . The bound in 
Eq. (|19|) allows for a variational determination of the trial 
propagator g(x,y), and the GEP is defined as the mini- 
mum of the functional 



VgbpM = V GE p[go} 



(21) 



with the optimal propagator go satisfying the gap equa- 
tion 



5V ( 



GEP 



Sg 



= o 



(22) 



f-Jo 



This stationary condition must hold at each point ip and 
the resulting propagator go depends on the value of ip. 
A trivial calculation yields 



V t 



GEP 



[g] = V c (ip) + V 1L (<p) + V [g] + V 2 [g] + V F [g}. (23) 



The first two terms are the classical potential and the 
one-loop term, and do not depend on the trial propagator 
since 



V 1L {<p) = N f (TrlogG m ) m=!w . 
The third term has no explicit dependence on ip 



V [g}=h\g} + ^-I${g}- 



d 4 P 



g(p) 

9b(p) 



(24) 



(25) 



Here g{jp) is the Fourier transform of g(x,y), the bare 
propagator is g^ 1 
defined as 



m 2 B + p 2 and the integrals Iq , 1\ are 

™=-\S (26) 



lo[g] = f-^g(p). 



(27) 



All the integrals in the four-dimensional Euclidean space 
are made finite by insertion of a cut-off A, and taking 
p 2 < A 2 . The fourth term V^g] has a quadratic explicit 
dependence on ip 



V 2 [g] 



Hg\ 



(28) 



and the last term Vf[<?] follows from the fermionic loop 
in the average of the quadratic part of Eq. (fT3"|) 



(29) 



where the kernel K(p) is the one- loop fermion polariza- 
tion function 



K(p) 



d 4 q 

W) 4 



Tr{G m (p + q)G m (q)} r 



(30) 



and G m (p) follows by Fourier transform of Eq. © 



G m (p) = (-TmP" + m)- 1 = 



p z + m z 



(31) 



Exact and approximate expressions for the kernel K(p) 
are reported in appendix A. 

According to Eq. ([22|) . the functional derivative of 
Eq. (|23"|) gives a gap equation that reads 



g^ =p 2 + n 2 [g a ] + N f y 2 K(p) (32) 
where the mass functional £l 2 [g] is defined as 



1 



1 



V M = m 2 B + -\ B ip 2 + -X B I [g] 



(33) 



and does not depend on p, while K{p) does not depend 
on g. In the limit y — > the optimal trial propagator 
takes the simple form g^ 1 = g^ 1 = p 2 + il 2 where the 
mass parameter Vt 2 = S! 2 [go] is the self-consistent solu- 
tion of Eq. (f3"3l . In this limit there is no wave function 
renormalization and we recover the GEP for a scalar the- 
ory 



V GE p(<fi) = V c {(p) + V a [g n ] + V 2 [gn] 



(34) 



where the last two terms have an implicit dependence on 
ip through gsi. This case has been studied in Ref.[12j in 
some detail. 

In general the propagator go has a more complicate 
dependence on p as a consequence of the Yukawa inter- 
action which adds the last term in Eq. (|32[) . We can study 
this dependence in some detail in the limit p — > 0. In this 
limit the kernel K(p) is an analytic function of p 2 and can 
be expanded in powers yielding 



K{p) « cto(m) + ai(m)p 2 + 0{p 4 ) 



(35) 



where the functions aj(m) are derived in Appendix A and 
must be evaluated for m = yip. The optimal propagator 
in Eq. (|32p can then be written in the same limit as 



9o(p) 



p 2 + M 2 R 



(36) 



where the wave function renormalization constant reads 

Z- 1 = 1 + y 2 N f a 1 {m) (37) 

and the renormalized mass is 

M 2 = Z [n 2 + y 2 N f a (m)] . (38) 

It is instructive to explore the content of the approx- 
imation in terms of graphs. The solution go of the gap 
equation Eq. (|32p can be regarded as the solution of the 
Dyson equation 



g(p) = gnip) + g(p) ■ [-y 2 N f K(p)\ ■ g n (p) 



(39) 



which sums up the whole class of ring-diagrams. From 
a formal point of view this kind of approximation is 
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equivalent to RPA, but here the self-consistent param- 
eter f2 2 = ri 2 [go] and the corresponding propagator 
g^ 1 = p 2 + fl 2 are functions of ip while in RPA they 
are fixed at their mean-field value. Thus the dependence 
of the effective potential on ip is different, as are the pre- 
dictions for the vev. 

Inserting back Eq. ([3"2"]) in Eq. (f2"3")) the effective potential 
takes the simple form 



V GE p(<p) = V c (<p) + V 1L {(p) + 



h - ^bIq 



(40) 



where we recognize the classic term, the one-loop correc- 
tion and the standard GEP term of a scalar 

theory [H 03 

(in brackets). However here the integrals Iq, 2j are eval- 
uated for g — go which is solution of the interacting gap 
equation Eq. (l32p . thus the last term does depend on the 
coupling y, and gives rise to large deviations from the 
simple one-loop result in the strong-coupling limit. 

The numerical evaluation of the effective potential can 
be obtained more easily by a simple integration of the 
first derivative 



dVoEP 9Vgep 



dp 



dip 



5V, 



GEP 



§9 



dp 



d 4 p 

WY 4 



(41) 



which by insertion of Eq. (|22|) is equivalent to the partial 
derivative 



dVoEP^p) dV G Ep{<p) 



dp 



dpi 



(42) 



yielding 
dV G 



EP 



dp 



m 2 B <p- 



3! V 



dVu 



dp 



-I ip- 



+ 



y 2 N f f d 4 p dK(p) 
(2tt) 4 dip 



g{p) 



(43) 



Here the trial propagator must be set to its optimal value 
go for any ip point, yielding an implicit dependence on 
that variable. Explicit expressions for the derivatives of 
K{p) and V\l are reported in appendix B and C respec- 
tively. 



III. COMPARISON WITH LATTICE DATA 

For any set of lagrangian parameters the GEP can be 
obtained by solving the gap equation Eq.(|32p and inte- 
grating the derivative in Eq. (|43p . In order to compare 
with lattice data some care must be taken in the choice 
of the energy units. While we would prefer units of the 
cut-off A, lattice data are usually reported in units of the 
lattice spacing a. Actually the cut-off can be thought as 
defining an effective lattice spacing a — c/A, with an un- 
known constant scale factor c = aA of order unity that 
depends on the approximation scheme and can be de- 
termined by a direct comparison with the lattice data. 




0.2 0.4 0.6 0.8 1 



m 



B 



Figure 1: The vev as a function of the bare mass ma in units 
of A/c for As = 0.1, y = 0.5, and for a scale factor c = 2.34 
(dashed line). With that choice of scale the GEP interpolates 
the lattice data of Ref.f2~Q| (circles). 



Once the constant scale factor is fixed, the predictions 
of the GEP can be compared with the lattice data(l2l|. 
Some lattice data on the Higgs-Top model are discussed 
in Ref.[21]. For a comparison we use the same set of 



parameters: y — 0.5, \b = 0.1 and 



= 0.1. As 



other lattice calculations, they take a large even number 
of fermions because of a well known duplication problem. 
We use the same value Nf = 8 in this section. The scale 
factor has been determined as c = 2.34 by a fit of the 
vev as shown in Fig.l. Here the vev is defined as the 
minimum point of the GEP, and is reported in units of 
A/c in order to fit the lattice data. 

With the scale factor fixed, the derivative of the ef- 
fective potential is evaluated by Eq. (|4"3")) and reported in 
Fig. 2 for the parameter set Xb = 0.1, y = 0.5, m 2 B = 0.1, 
Nf = 8, in units of the effective lattice spacing a = c/A. 
The lattice data of Ref . j2l[ are shown in the same figure 
for comparison. We find a good general agreement, with 
no need for any other tuning of the parameters. Unfortu- 
nately we only found detailed lattice data for comparison 
in this weak-coupling limit where the variational method 
does not add to much to the usual perturbative treat- 
ment of the model. Some recent lattice data have been 
reported in the strong coupling limit but for the special 
limit cases Xb — > oo and Xb = o[24[. Strong-coupling de- 
viations from the one-loop predictions will be discussed 
in the next section. 

However, even in this weak coupling limit, at variance 
with the perturbative method that predicts an unsta- 
ble vacuum for large values of ip, the GEP is perfectly 
bounded and predicts a stable vacuum as shown in Fig. 3 
in agreement with the lattice calculations, confirming 
that the instability of the vacuum is not a consequence 
of the interaction with the fermions, but a misleading 
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0.15 r 
0.1 - 




-0.15 1 1 1 1 1 1 

0.5 1 1.5 2 

9 

Figure 2: The derivative of the effective potential dV/dtp is 
reported as a function of the variable tp in units of A/c with 
c = 2.34, and for the parameter set As = 0.1, y = 0.5, m% = 
0.1, Nf = 8. The GEP (solid line) is compared with the 
lattice data points of Ref.[2l| (error bars). 




2 4 6 8 10 12 

9 

Figure 3: The derivative of the effective potential dV/dtp is 
reported as a function of the variable tp in units of A/c with 
c = 2.34, and for the same parameter set as for Fig. 2, but for 
a wide range of tp. The GEP (solid line) is compared with 
the lattice data points of Ref.|2ll| (circles). The derivative 
is positive for large values of tp, even when tp > A, and no 
vacuum instability occurs in the GEP 

internal problem of the standard perturbative approxi- 
mation, emerging when the field is much larger than the 
cut-off|21]. Actually vacuum instability emerges in the 
renormalized perturbation theory because the renormal- 
ized couplings are allowed to assume any value, whereas 
in bare theories, like the present variational method or 
lattice calculations, only a limited range of renormalized 
couplings are possible, which vanish logarithmically with 




-0.34 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 
0.5 1 1.5 2 2.5 3 3.5 4 

9 

Figure 4: The GEP VGEp(tp) is shown as a function of the 
variable tp in units of the effective lattice spacing a — c/A 
with c = 2.34, and for the same parameter set as for Fig. 2. 

the cut-off. By inspection of Eq. (|4"3"|) we observe that in 
the unphysical limit of tp 3> A the integrals become irrel- 
evant as they are cut at a relatively small value of A. The 
negative unstabilizing one- loop term vanishes as ~ 1/tp 
in that limit according to Eq. (|C3[) , and the leading term 
in Eq. (|43p is the positive classic tp 3 term. Thus the ef- 
fective potential cannot be unbounded for large values of 
the field. 

For the same set of parameters the GEP is reported in 
Fig. 4 in units of the effective lattice spacing. The min- 
imum is at tp = v w 2.1 and in the same units the top 
mass is m — yv ~ 1. It is remarkable that the breaking 
of symmetry is entirely driven by the Yukawa interaction 
with the fermions: here m B > 0, the classical potential 
is symmetric and would predict a vanishing vev. Thus 
in the SM the symmetry is broken because of the large 
Yukawa coupling of the top quark y « 0.7 that gives rise 
to a finite vev and gives back a mass to all the fermions. 
There is no need for unphysical negative bare squared 
masses m? B < in the Lagrangian, that are quite hard 
to be explained for a free field theory. Moreover in sev- 
eral extensions of the SM the Higgs sector can be quite 
simplified, because there is no need to assume that the 
classical potential has a mimimum at a broken symmetry 
point ip ^ 0. For instance, the minimal left-right sym- 
metric extension of the SM(25l [26| does not require the 
existence of scalar bi-doublet fields, and a model with 
only two doublets is perfectly viable [13, [Hj]. 



IV. STRONG COUPLING LIMIT 

The GEP has already been studied for large values of 
the self-coupling of the scalar field in the past [HI]. 
The method has been shown to be reliable for large cou- 
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plings, even if these are almost ruled out by the recent 
experimental evidence of a light Higgs mass. In this sec- 
tion we would like to explore the limit of a strong Yukawa 
coupling, much larger than Xb- We take = 0.1, and 
Nf = 3 in order to represent a QCD quark multiplet like 
the top quark. Moreover we use units of A in this section 
(c = 1), and explore the model for a large coupling y and 
generic values of the bare mass. Thus in these units we 
are left with two free parameters, y and m 2 B . However, 
the phyisical requirement of a non vanishing small top 
mass, quite smaller than the cut-off (m <C 1 in our units), 
decreases our degree of freedom, and limits the ratio of 
the free parameters in the range 0.2 < ms/y < 0.25. In 
fact the physical masses can be regarded as the inverse 
of correlation lengths £ = 1/m, and the condition m <C 1 
is equivalent to the requirement that £ is much larger 
than the effective lattice spacing. Only in that limit the 
model makes sense. Actually, we can show that the pre- 
dictions of the one-loop approximation only depend on 
the single parameter y/m,B in the strong coupling limit 
!/ > Ab » 1. When the physical top mass m = yv is dis- 
played as a function of ms/y, the one-loop data collapse 
on a single curve in the strong coupling limit, and do not 
depend on y as dispayed in Fig. 5. 



tential follows from Eq.(|43p 
1 dV G EP 



dm 



\ dm 2 



Xj 



X B 
¥ 7 ° 



y 2 N f 



m B 
V 

d 4 p dK(p) 
(2tt) 4 dm 2 



B 3 

3!y 4 



S(p). (44) 



We recognize a first term in square brackets arising from 
the classical plus one-loop potential; a second term of or- 
der ~ Xs/y, arising from the standard GEP correction 
to the scalar potential; a third term of order ~ y 2 , aris- 
ing from our non-perturbative treatment of the Yukawa 
interaction. We observe that V\l and K have no explicit 
dependence on y and ms when written as functions of 
the variable m. Thus at one-loop, in the strong coupling 
limit y Xb, the derivative of the effective potential 
reads 



dV_ 

dm 



\ dm 2 



m B 

y 



(45) 



and only depends on the ratio ms/y- The stationary 
points of the potential are the solution of the equation 
dV/dm = 0. When the symmetry is broken we find a 
maximum at m — and a minimum given by the solution 
of the equation 
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Figure 5: The top mass as a function of the parameter 
(ms/y), for a range of Yukawa couplings going from y = 3 
to y = 20, according to the GEP. Energies are in units of 
the cut-off A. The one-loop data collapse on a single curve 
(dashed). 
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(46) 



thus a plot of the one-loop top mass m at the minimum, 
as a function of the single parameter (mB/y), must follow 
a single curve for any strong coupling y. 
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Evidence for this strong coupling scaling can be found 
by a simple analysis of the single terms contributing to 
the effective potential. By taking the mass m = yip as 
independent variable, the derivative of the effective po- 



Figure 6: The top mass as a function of the Yukawa coupling 
y, for a typical value of the parameter (ms/y) = 0.21. The 
prediction of the GEP (solid circles) is compared with the 
steady one- loop data (open circles). Masses are in units of 
the cut-off. 
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Figure 7: The wave function renormalization constant Z is 
reported according to Eq. (|37[) as a function of the Yukawa 
coupling, for a typical value of the parameter (ms/j) = 0.21. 



This single curve predicts a universal critical point 
(ms/y) c for the single parameter, when the maximum 
and minimum coincide: 



\dm 2 /B ,_„ 



1 ( HL£ 

2 >v y 



(47) 



and from the explicit expression in Eq. (|C3|) we obtain 



m B 

y 



2tt 



(48) 



that for Nf — 3 is (ms/j/) c ~ 0.276. In this strong cou- 
pling regime, deviations from the single parameter scal- 
ing can only arise from the non-perturbative last term in 
Eq. (|44|l . Thus deviations from scaling can be used as a 
measure of the non-perturbative effects in the GEP. The 
top mass m, as emerging from the GEP calculation, is 
reported in Fig. 5 as a function of the scaling parameter 
itib /y for several values of the Yukawa coupling ranging 
from y = 3 to y = 20. All the one-loop data collapse on 
the dashed line. Deviations from the scaling are present 
when the GEP is considered, but they are small for the 
phenomenological value y 0.7. These deviations can- 
not be neglected in the range of large couplings. In Fig. 6 
the top mass is reported as a function of the Yukawa cou- 
pling y for a typical value of the parameter ms/y = 0.21. 
We notice the different behaviour of one-loop and GEP 
approximations: the top mass is reduced by fluctuations. 
Moreover the fluctuations give rise to a decreasing of the 
wave function renormalization constant Z according to 
Eq. (|37p as shown in Fig. 7. This effect is negligible for 
weak couplings, but becomes very large in the strong 
coupling limit. 

More insights come from a closer study of the Higgs 
propagator go(p) that solves Eg. (f3"2"|) . The renormalized 



Figure 8: The renormalized inverse propagator Zg~ x is shown 
as a function of p 2 , in units of the cut-off A, for a typical 
value of the parameter (mj/j) = 0.21, and for a moderate 
coupling y — 3 (solid line). Here g is the optimal propagator 
solving Eq. (|32p . The linear approximate inverse propagator 
of Eq. (|36[) is shown for comparison (dashed line). 
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Figure 9: The same as Fig. 8 but for a strong Yukawa coupling 
y = 20 (solid line). The linear approximate inverse propagator 
is shown for comparison (dashed line). 



inverse propagator Zg^ 1 is reported as function of p 2 for 
a typical value of the parameter ms/y — 0.21 in Figs. 8 
and 9. The approximate linear expression in Eq. (|36[) . 
derived for small values of p, is shown for comparison. 
In Fig. 8, for a moderately large Yukawa coupling y = 3, 
deviations from linearity are small and only occur for 
very large values of p 2 . Deviations are large in the strong 
coupling limit, as shown in Fig. 9 for y = 20. 

A shortcoming of the approximation comes from the 
break-down of the convergence criterion when approach- 
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ing the critical point. Before reaching the critical point 
the fluctuations of the Higgs field become large enough 
to make the expansion in Eq.(fT0|) useless. In fact in Fig. 5 
the plots do not reach the critical point, but we checked 
that the criterion of convergence is satisfied for the re- 
ported data. From a technical point of view the renor- 
malized mass in Eq. (|38l) becomes small when approach- 
ing the critical point, and it vanishes before reaching the 
transition. A real pole occurs in the Higgs propagator 
and the integrals diverge. However this is just a sign that 
the Higgs fluctuations are too large and the expansion in 
Eq. pop does not hold any more. In our calculation we 
never let the renormalized mass be too small, in order to 
fullfill the convergence criterion. Thus our data are still 
reliable even if they cannot span the entire range of the 
free parameters. 



V. CONCLUSIONS 

We have studied the Higgs- Top model by a non- 
perturbative variational extension of the GEP that in- 
corporates fermions. While the pure GEP is known to 
give trivial results for the fermions[2(|, some effects of 
their fluctuations have been included in the GEP by an 
hybrid method: fermions are integarted out exactly and 
the resulting effective action is expanded in powers of the 
fluctuating Higgs field. Even in the strong coupling limit 
a second order expansion provides reliable predictions for 
a large range of the free parameters. By fixing an effec- 
tive lattice spacing, the GEP is found in good general 
agreement with the available lattice data for the model. 

In the strong coupling limit the gaussian fluctuations 
reduce the vev and the top mass, as displayed by a com- 
parison with the one-loop approximation. In this limit 
the one-loop data are seen to follow a single-parameter 
scaling, while the gaussian fluctuations give rise to im- 
portant deviations from scaling. 

At variance with renormalized perturbation theory, 
this study confirms that the Higgs- Top model has a per- 
fectly stable vacuum even in the strong coupling limit 
and that, in presence of Yukawa couplings with fermions, 
there is no need for a finite vev in the classical potential. 
In fact a non-vanishing vev can be predicted even when 
m 2 B > and the classical potential has no broken sym- 
metry vacuum. Thus fluctuations are very relevant for 
determining the correct vev, and naive discussions based 
on the classical potential cannot be used for the Higgs- 
Top model. 

On general grounds, inclusion of fermions in a non- 
perturbative and almost analytical calculation is of inter- 
est by itself. Moreover the same method could be used in 
the framework of effective theories [29j or in a low-energy 
study of gauge theories like QCD when the perturbative 
approximation breaks down. While the GEP was derived 
for the SU(2) gauge theory in the past(l6j. no previous 
attempt has been reported for incorporating fermions. 
Thus the present study could be extended in that direc- 



tion. 



Appendix A: Polarization function K(p) 

The fermionic polarization function Kip) has been 
calculated exactly, integrating inside the hyper-sphere 
p 2 < A 2 in the four-dimensional Euclidean space. In- 
serting Eq.(|3ip in Eq. (13011 and evaluating the trace 



K(p) 



m 2 — q 2 — q ■ p 



l A (27r) 4 (m 2 + q 2 )(m 2 + q 2 + p 2 + 2q ■ p) ' 

(Al) 

Here the Feynman trick cannot be used because of the 
finite cut-off. However by a tedious but straightforward 
calculation the exact polarization function can be written 
as 

K{p) = K + K^p) + K 2 {p) + K 3 (p) 
where Kq is the constant 
1 



K = — 



8tt 2 



a 2 2 , A + ™> 

A — m log - — 



(A2) 



(A3) 



and the three functions Ki(p) follow 
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K 3 (p) 



4m 2 +p 2 T , . 
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where ft and 8 are the functions 

(3(p) — m 2 + A 2 — p 2 



S(p) = \/l3 2 + Am^p 



2„2 



and J is the integral 



J(P) 



dx 



yjx — p 2 ) 2 + im 2 p 2 



(A4) 
(A5) 

(A6) 

(A7) 
(A8) 

(A9) 



The integral J has been evaluated in Ref.[12] and its 
explicit expression is 



r, \ /3 + 6 -2m 2 2 
J(P) = o p 9 

where p is the function 

to / 1 1 \ 



(A10) 



Pip) = log 



p \t 2 ti 



Am 2 



(t 2 - *-)(tl - *+) 



(ta -*+)(*! -*-) 



(All) 
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and its arguments are 



2m 



±\ 1 



4m 2 



(A12) 



riif 3 (p) = (4m 2 + p 2 ) dJ(p) _ J(p) 
dm 2 167r 2 p 2 dm 2 47r 2 p 2 ' 



(B6) 



h = 



t 2 = 



p 



ft + 5 
2mp 



(A13) 



It can be easily seen that the singular terms cancel in 
Eq. (|A2[) and the resulting function K(p) is an analytic 
function of p 2 . Its expansion in powers of p 2 reads 



K{p) 



n=0 



a n (m)p 2 



(A14) 



where the coefficients a n (m) are functions of m 2 . In units 
of the cut-off A the first two coefficients of the expansion 
are 



a (m) 



1 



1 



2m 2 



1 + m 2 



3m log 



1 + m 2 



(A15) 



where the functions ft, 5 and J are defined in Appendix 
A. The derivative of J follows as 



dJ(p) ft + 2p 2 -5 2p 2 (A 2 + 2m 2 ) 
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. 1 1 + m 2 (6m 4 + 21m 2 + 7) 
a l( m) = — log — 48 ,2( 1+m2) 3 ■ ( Al6 ) 



where the functions t±, t\, t<z are defined in Appendix A. 



Appendix B: Derivative of K(p) 

Explicit expressions for the exact and approximate po- 
larization function K(p) have been reported in Appendix 
A. The derivative can be written as 



where 

dK(p) 
dm 2 
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and the single terms Ki are given in Appendix A. By an 
explicit calculation we obtain 
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Appendix C: One-loop potential V\l and its 
derivative 



The one- loop term Viz, (92) is defined in Eq. (j2~4")) and 
its explicit evaluation, in terms of m = yip, yields the 
integral 

V 1L [m) = -2N f £ 2 -^^\og(m 2 + p 2 ) (CI) 

and in units of the cut-off A 
Nt 
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The derivative follows 
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